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The recent detection of blazar 3C279 by MAGIC has confirmed previous indications by 
H.E.S.S. that the Universe is more transparent to very-high-energy gamma rays than pre- 
viously thought. We show that this fact can be reconciled with standard blazar emission 
models provided photon oscillations into a very light Axion-Like Particle occur in ex- 
tragalactic magnetic fields. A quantitative estimate of this effect explains the observed 
spectrum of 3C279. Our prediction can be tested in the near future by the satellite-borne 
GLAST detector as well as by the ground-based Imaging Atmospheric Cherenkov Tele- 
scopes H.E.S.S., MAGIC, CANGAROO III, VERITAS and by the Extensive Air Shower 
arrays ARGO-YBJ and MILAGRO. 

1 Introduction 

As is well known, in the very-high-energy (VHE) band above 100 GeV the horizon of the 
observable Universe rapidly shrinks as the energy further increases. This comes about because 
photons from distant sources scatter off background photons permeating the Universe, thereby 
disappearing into electron-positron pairs [T]. The corresponding cross section (7(77 — > e + e~) 
peaks where the VHE photon energy E and the background photon energy e are related by e ~ 
(500 GeV/E) eV. Therefore, for observations performed by Imaging Atmospheric Cherenkov 
Telescopes (IACTs) - which probe the energy interval 100 GeV — 100 TeV - the resulting cosmic 
opacity is dominated by the interaction with ultraviolet/optical/infrarcd diffuse background 
photons (frequency band 1.2 • 10 3 GHz — 1.2 • 10 6 GHz, corresponding to the wavelength range 
0.25 /Ltm — 250 /jm), usually called Extragalactic Background Light (EBL), which is produced by 
galaxies during the whole history of the Universe. Neglecting evolutionary effects for simplicity, 
photon propagation is controlled by the photon mean free path \y(E) for 77 — > e + e~, and so 
the observed photon spectrum <J> b s (-E, D) is related to the emitted one & em (E) by 

<S> oha {E,D)=e- D ' x ^ $ cm (£) . (1) 

Within the energy range in question, \ 7 (E) decreases like a power law from the Hubble 
radius 4.2 Gpc around 100 GeV to 1 Mpc around 100 TcV %. Thus, Eq. (TJ entails that the 
observed flux is exponentially suppressed both at high energy and at large distances, so that suf- 
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ficiently far-away sources become hardly visible in the VHE range and their observed spectrum 
should anyway be much steeper than the emitted one. 

Yet, observations have not detected the behaviour predicted by Eq. (p}. A first indication 
in this direction was reported by the H.E.S.S. collaboration in connection with the discovery 
of the two blazars H2356-309 (z = 0.165) and 1ES1101-232 (z = 0.186) at E ~ 1 TeV [3]. 
Stronger evidence comes from the observation of blazar 3C279 (z = 0.536) at E ~ 0.5 TeV by 
the MAGIC collaboration g]. In particular, the signal from 3C279 collected by MAGIC in the 
region E < 220 GeV has more or less the same statistical significance as the one in the range 
220 GeV < E < 600 GeV {6.1a in the former case, 5.1a in the latter). 

A suggested way out of this difficulty relies upon the modification of the standard Synchro- 
Sclf-Compton (SSC) emission mechanism. One option invokes strong rclativistic shocks [5]. 
Another rests upon photon absorption inside the blazar [6]. While successful at substantially 
hardening the emission spectrum, these attempts fail to explain why only for the most distant 
blazars does such a drastic departure from the SSC emission spectrum show up. 

Our proposal - usually referred to as the DARMA scenario - is quite different [7] . Implicit in 
previous considerations is the hypothesis that photons propagate in the standard way through- 
out cosmological distances. We suppose instead that photons can oscillate into a new very light 
spin-zero particle - named Axion-Likc Parlicle (ALP) - and vice- versa in the presence of cosmic 
magnetic fields, whose existence has definitely been proved by AUGER observations [5]. Once 
ALPs are produced close enough to the source, they travel unimpeded throughout the Universe 
and can convert back to photons before reaching the Earth. Since ALPs do not undergo EBL 
absorption, the effective photon mean free path A 7i eff(-E) gets increased so that the observed 
photons cross a distance in excess of A 7 (-E). Correspondingly, Eq. {T]) becomes 

<& obs (£, D) = e- fl A„**(*) $ em (E) , (2) 

from which we see that even a slight increase of A 7i eff(-E) gives rise to a huge enhancement of 
the observed flux. It turns out that the DARMA mechanism makes A 7je ff(i') shallower than 
A 7 (i?) although it remains a decreasing function of E. So, the resulting observed spectrum is 
much harder than the one predicted by Eq. |T]) , thereby ensuring agreement with observations 
even for a standard SSC emission spectrum. As a bonus, we get a natural explanation for the 
fact that only the most distant blazars would demand <f> cm (i?) to substantially depart from the 
emission spectrum predicted by the SSC mechanism. 

Our aim is to review the main features of our proposal as well as its application to blazar 
3C279. 

2 DARMA scenario 

Phenomenological as well as conceptual arguments lead to view the Standard Model of particle 
physics as the low-energy manifestation of some more fundamental and richer theory of all 
elementary-particle interactions including gravity. Therefore, the lagrangian of the Standard 
Model is expected to be modified by small terms describing interactions among known and new 
particles. Many extensions of the Standard Model which have attracted considerable interest 
over the last few years indeed predict the existence of ALPs. They are spin-zero light bosons 
defined by the low-energy effective lagrangian 

Calf = \d» ad^a-^-m 2 a 2 - -J-F F^a , (3) 
A A 
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where F^ v is the electromagnetic field strength. is its dual, a denotes the ALP field whereas 

— 1/2 

m stands for the ALP mass. According to the above view, it is assumed M 3> G F ~ 250 GeV. 
On the other hand, it is supposed that m <C G F ~ 250 GeV. The standard Axion |Hj is the 
most well known example of ALP. As far as generic ALPs are concerned, the parameters M 
and m are to be regarded as independent. 

So, what really characterizes ALPs is the trilincar 7-7-a vertex described by the last term 
in £alP; whereby one ALP couples to two photons. Owing to this vertex, ALPs can be emit- 
ted by astronomical objects of various kinds, and the present situation can be summarized as 
follows. The negative result of the CAST experiment designed to detect ALPs emitted by the 
Sun yields the bound M > 0.86 ■ 10 10 GeV for m < 0.02 eV [10]. Moreover, theoretical consid- 
erations concerning star cooling via ALP emission provide the generic bound M > 10 10 GeV, 
which for m < 10~ 10 eV gets replaced by the stronger one M > 10 11 GeV even if with a large 
uncertainty The same 7-7-a vertex produces an off-diagonal element in the mass matrix 
for the photon- ALP system in the presence of an external magnetic field B. Therefore, the in- 
teraction eigenstatcs differ from the propagation eigenstates and photon- ALP oscillations show 

up ng. 

We imagine that a sizeable fraction of photons emitted by a blazar soon convert into ALPs. 
They propagate unaffected by the EBL and we suppose that before reaching the Earth a sub- 
stantial fraction of ALPs is back converted into photons. We further assume that this photon- 
ALP oscillation process is triggered by cosmic magnetic fields (CMFs), whose existence has 
been demonstrated very recently by AUGER observations [8j. Owing to the notorious lack 
of information about their morphology, one usually supposes that CMFs have a domain-like 
structure |13| . That is, B ought to be constant over a domain of size Ldom equal to its coher- 
ence length, with B randomly changing its direction from one domain to another but keeping 
approximately the same strength. As explained elsewhere [H] , it looks plausible to assume the 
coherence length in the range 1 — lOMpc. Correspondingly, the inferred strength lies in the 
range 0.3 - 1.0 nG [Hj. 

3 Predicted energy spectrum 

Our ultimate goal consists in the evaluation of the probability P 1 ^ 1 {E,D) that a photon 
remains a photon after propagation from the source to us when allowance is made for photon- 
ALP oscillations as well as for photon absorption from the EBL. As a consequence, Eq. (0) 
gets replaced by 

$ ob8 (E, D) = P 7 ^ 7 (E, D) $ em (E) . (4) 

We proceed as follows. We first solve exactly the beam propagation equation arising from £alp 
over a single domain, assuming that the EBL is described by the "best-fit model" of Knciske 
et al. [15) . Starting with an unpolarized photon beam, we next propagate it by iterating the 
single-domain solution as many times as the number of domains crossed by the beam, taking 
each time a random value for the angle between B and a fixed overall fiducial direction. We 
repeat such a procedure 10 000 times and finally we average over all these realizations of the 
propagation process. 

We find that about 13% of the photons arrive to the Earth for E = 500 GeV, representing 
an enhancement by a factor of about 20 with respect to the expected flux without DARMA 
mechanism (the comparison is made with the above "best-fit model"). The same calculation 
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gives a fraction of 76% for E = 100 GeV (to be compared to 67% without DARMA mechanism) 
and a fraction of 3.4% for E = 1 TeV (to be compared to 0.0045% without DARMA mechanism) . 
The resulting spectrum is exhibited in Fig. 1. The solid line represents the prediction of the 
DARMA scenario for B ~ 1 nG and idom — 1 Mpc and the gray band is the envelope of the 
results obtained by independently varying B and L& om within a factor of 10 about such values. 
These conclusions hold for m <C 10 -10 eV and we have taken for definitcness M ~ 4 ■ 10 11 GeV 
but we have choked that practically nothing changes for 10 11 GeV < M < 10 13 GeV. 

Our prediction can be tested in the near future by the satellite-borne GLAST detector as 
well as by the ground-based IACTs H.E.S.S., MAGIC, CANGAROO III, VERITAS and by the 
Extensive Air Shower arrays ARGO-YBJ and MILAGRO. 




Photon energy (GeV) 



Figure 1: The two lowest lines give the fraction of photons surviving from 3C279 without the 
DARMA mechanism within the "best-fit model" of EBL (dashed line) and for the minimum 
EBL density compatible with cosmology (dashed-dotted line) [15] . The solid line represents the 
prediction of the DARMA mechanism as explained in the text. 
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